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1 INTRODUCTION 



Recently, [Avenhaus, Schmid & Meyer| ( |2012[ > conducted a 
search for debris discs in a sample of 85 M stars using data from 
the Wide-field Infrared Survey Explorer (WISE). They reported a 
null detection rate, which they attributed to either the different evo- 
lution of dust around M stars or to an age effect. By contrast, debris 
disc detections are prevalent around A, F and G stars ( |Rieke et"ar 
2005} |Bry den etal. (2 006 Chen et al.|2006l|Su et al.|2006||Wyatt 



et al.|2007 



Meyer et al. 2008, Urb an et al.|20I2| >. The null result 
of |Avenhaus, Schrmd & Meyer ( 2012) motivates a theoretical inter- 
pretation in terms of the plcmetesimah that are widely believed to 
reside in these discs (Krivov et al. 2008, Hens & Tremaine 2010, 
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ABSTRACT 

Motivated by the reported dearth of debris discs around M stars, we use survival models to 
study the occurrence of planetesimal discs around them. These survival models describe a 
planetesimal disc with a small number of parameters, determine if it may survive a series of 
dynamical processes and compute the associated infrared excess. For the WISE satellite, we 
demonstrate that the dearth of debris discs around M stars may be attributed to the small semi- 
major axes generally probed if either: 1 . the dust grains behave like blackbodies emitting at a 
peak wavelength coincident with the observed one; 2. or the grains are hotter than predicted by 
their blackbody temperatures and emit at peak wavelengths that are shorter than the observed 
one. At these small distances from the M star, planetesimals are unlikely to survive or persist 
for time scales of 300 Myr or longer if the disc is too massive. Conversely, our survival models 
allow for the existence of a large population of low-mass debris discs that are too faint to be 
detected with current instruments. We gain further confidence in our interpretation by demon- 
strating the ability to compute infrared excesses for Sun-like stars that are broadly consistent 
with reported values in the literature. However, our interpretation becomes less clear and large 
infrared excesses are allowed if only one of these scenarios holds: 3. the dust grains are hotter 
than blackbody and predominantly emit at the observed wavelength; 4. or are blackbody in 
nature and emit at peak wavelengths longer than the observed one. Both scenarios imply that 
the parent planetesimals reside at larger distances from the star than inferred if the dust grains 
behaved like blackbodies. In all scenarios, we show that the infrared excesses detected at 22 
/xm (via WISE) and 70 jiva (via Spitzer) from AU Mic are easily reconciled with its young age 
(12 Myr). Conversely, the existence of the old debris disc (2-8 Gyr) from GJ 581 is due to the 
large semi-major axes probed by the Herschel PACS instrument. We elucidate the conditions 
under which stellar wind drag may be neglected when considering dust populations around 
M stars. The WISE satellite should be capable of detecting debris discs around young M stars 
with ages ^10 Myr 

Key words: minor planets, asteroids - planets and satellites: general 



The traditional view of debris discs is that they are dust discs 
generated by collisions between parent planetesimals, which are 
much larger in size ( [Backman & Paresce|1993| |Zuckerman|2001[ 
|Wyatt|[2008^ . Only dynamically hot discs produce debris discs. 
These are discs in which collisions occur and produce enough 
dust such that the infrared emission associated with the dust grains 
(from reprocessing the incident starlight) exceeds the detection 
threshold and is hence observable using an infrared telescope. In 
other words, debris discs are the progeny of dynamically hot plan- 
etesimal discs. By contrast, dynamically warm discs contain plan- 
etesimals with smaller sizes and in larger numbers than those found 
in dynamically hot discs ( |Heng & Tremaine|2010| ). Collisions occur 
frequently between the planetesimals, but they are non-destructive 
and do not produce dust in any significant amount. However, the 
planetesimals are small and numerous enough that they provide a 
non-negligible covering fraction around the star. It is the planetes- 
imals themselves that reprocess the incident starlight into infrared 
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Figure 1. Distances from the star probed as a function of tlie stellar type 
for the four WISE wavebands, estimated using equation j21) . For compar- 
ison, we include the distances associated with the Spitzer 70 fj,m channel. 
Note that we have assumed /dust = 1 in these plots. For small dust grains 
(/dust > 1). the distances involved are larger and the chances of survival 
are more favourable; in this context, these distances are lower limits. See 
§3.2|for a discussion of other interpretations when /dust ¥" 1- 



flux. In essence, dynamically warm discs may mimic debris discs, 
although most of the debris discs observed so far are not believed 
to be dynamically warm discsrj 

There are two plausible approaches to modeling dynamically 
hot and warm discs. The first approach is to construct a formation 
model, starting from the birth of the gaseous, dusty, protoplanetary 
disc and ending with a gas-poor disc populated with planetesimals 
(and possibly planets). The multitude of poorly-understood pro- 
cesses associated with planet formation renders this a formidable 
task. Parametrizing our ignorance of these processes produces a 
model with a large number of free parameters. The second ap- 
proach is to construct a survival model jHeng & Tremaine|2010| ), 
which addresses the following question: given a planetesimal disc 
described by a small number of parameters, can it survive a series 
of dynamical processes that act to destroy it on time scales possibly 
smaller than its age? 

In the present study, we employ the survival models of |Heng| 
|& Tremaine| ( |20I0l l to examine a suite of model discs that are both 
dynamically hot and warm. We apply these models to discs around 
M stars and specialize to the wavebands of WISE. We demonstrate 
that a possible reason for the dearth of debris discs observed around 
M stars is because the WISE wavebands generally probe semi- 
major axes ~ 1 AU or less (Figure [TJ, where dynamical survival 
conditions are unfavourable for the parent planetesimals in discs to 
persist for ~ 300 Myr or longer. This interpretation becomes less 
clear when the associated dust grains or planetesimals deviate from 



^ One way of breaking the degeneracy in interpretation is to search for "sil- 
icate features" in the spectral energy distribution of the observed disc. For 
example, HD 72095, HD 69830 and r/ Corvi are known to possess such 
silicate features, which provides direct evidence for the presence of small 
dust grains (see |Heng & Tremaine|20f0, and references therein). Another 
method is to compare the disc size, inferred from fitting a blackbody to the 
excess emission, with the actual size if the disc is resolved; if the resolved 
size is larger, then small dust grains are likely to be dominating the emis- 
sion. Yet another test is to determine if the far-infrared and sub-millimetre 
spectrum of the disk follows a blackbody or exhibits a steeper drop-off due 
to inefficient emission from small grains, which argues for the presence of 
a dynamically hot disc. 



blackbody behaviour or emit at wavelengths longer than 22 fira. 
We gain further confidence in our interpretation by demonstrating 
that we can produce infrared excesses for discs around Sun-like 
stars that are broadly consistent with reported values ( Bryd en et al.| 
|2006[[Chen et al..2006j , as well as obtain model disc solutions for 
the infrared excesses observed from AU Mic and GI 581. 

In !j2] we review the observational constraints. In ip] we con- 
struct our survival models of both dynamically hot and warm plan- 
etesimal discs. We present our results in sQand discuss their impli- 
cations in Sp] 



2 OBSERVATIONAL CONSTRAINTS 



[Avenhaus, Schmid & Meyer|( |20I2| l searched for debris discs 
at 3.4 fjm and 12 pm around 85 M stars and at 22 ^m around 84 
M stars using data from the WISE satellite. They reported a de- 
tection rate of 0.01qq% and estimated the detection threshold for 
the infrared excess to be /la.thros ~ 0.1-0.3. We define the in- 
frared excess to be /ir = {Fx,obs — -Fa,*)/-Fa,* where FA.obs is 
the observed flux at a wavelength A and Fa,* is the photospheric 
flux from the star at the same wavelength (as given by, e.g., Kurucz 
models). While no stellar ages are tabulated, [Avenhaus, Schmid &| 
Meyer (2012) do remark that they expect the majority of the stars 
in their sample to have ages exceeding 300 Myr There is little age 
information and the sample is for nearby M stars. (See also ,Gautier| 
let al.|2008| and |Simon et al.|2012| ) 

Additionally, [Avenhaus, Schmid & Meyer|p012 l) reported a 
tentative infrared excess for the relatively bright AU Micr|(Ml star, 
t* ~ 12 Myr), part of the /3 Pic moving group, of /ir — 0.067 ± 
0.02 at 22 ^m. It was previously known that AU Mic has an infrared 
excess of /iR = 13.1 ± 1.5 at 70 ^m jPlavchan et al.[2009[ ). 

For comparison, we note that ~ 10% of Sun-like stars, with 
ages ~ 10 Myr to ~ 1 Gyr, have reported infrared excesses of 
/iR ~ l-IO ( [Chen et al.|"2 006; Bryden e t al.[[2006[ [Meyer et al.[ 
[2008[ l. For example, [Meyer et al., (|2008J find the debris disc fre- 
quency around Sun-like stars to be < 4% for A = 24 fira and 
t* > 300 Myr. Most systems with detected 24 /im excesses are in- 
ferred to peak at longer wavelengths. Generally, infrared excesses 
are more commonly detected at 70 fj,m. (As an aside, we note that 
[Rizzuto, Ireland & Zucker|2012| detected 22 /im infrared excesses 
from B, A and F stars in the 5-17 Myr-old Sco-Cen association 
with WISE, which they interpreted to emanate from debris discs.) 



3 METHOD 

Our model contains 9 basic parameters. 

(i) Star: the stellar mass Ali,, age f*, effective temperature T* 
and radius 7?* . 

(ii) Disc: the total mass i\fdisc and semi-major axis a of the disc. 

(iii) Planetesimals: the radius r, radial velocity dispersion a^ 
and mass density p of the planetesimal. 

The stellar luminosity is given by £* — A-KR^assT^ where ass 
is the Stefan-Boltzmann constant. We obtain the values of the stel- 
lar parameters from Table 1 of [Kaltenegger & Traub[ ( [2009[ ). We 
assume p = 3 g cm~^. 



Also known as HD 197481 and GJ 803. 
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3.1 Survival Models for Dynamically Hot & Warm Discs 

In this section, we state the conditions for the survival mod- 
els of both dynamically hot and warm discs. Many of these condi- 
tions are previously described in |Heng & Tremainel ( |2010| ), but we 
rewrite them in the form of inequalities involving r and/or Or and 
also generalize them for A/* 7^ Mq . 

We examine a disc that is centered at a distance a from its star 
and extends from a — fmCi/i to a + fmCi/i, i.e., it has a width of 
fma/2. When fm = 1, we term the disc to extend over an "octave" 
following the terminology of|Heng & Tremainel(|2010^. 



3.1.1 Common Conditions 

A trivial condition for a disc with a mass A/disc is that it needs 
to contain more than one planetesimal. 



r < 



\ 47rp 



(1) 



We demand that the planetesimal eccentricities and inclina- 
tions are not too large: eo < /e where eo is the root-mean-square 
(rms) eccentricity and fe = 0.5. We assume that io/co = 0.5 
where ig is the rms inclination. Requiring eg < fe is equivalent to 
demanding that the disc is thin. 



o-,. < fc 



V 2a 



1/2 



(2) 



A disc is considered to be dynamically hot or warm when the 
radial excursions of the planetesimals (2aeo) exceed the typical ra- 
dial separation between planetesimals, 



ar > 



7i-/mpr3 fGAU 



3Mdisc V 2a 



1/2 



(3) 



If the planetesimals are numerous enough that the disc may 
be approximated as a fluid, then the gravitational stability of a hot 
disc is described by the usual Toomre criterion ( |Toomre||1964| ). 
However, if the fluid approximation breaks down, then the criterion 
needs to be generalized as described in §3.2.1 of |Heng & Tremaine] 
pOTOf. Specifically, this occurs when 

rX^r'V^V'. (4) 



\pALj V /" 
When equation iHb is fulfilled, gravitational instability occurs when 



2 ^ 



Gpr^ 



pr-^M^, ( f, 



Md 



Otherwise, the usual Toomre criterion applies. 



A/dis. 



G 



1/2 



(5) 



(6) 



3.1.2 Dynamically Hot Discs 

A dynamically hot disc necessarily needs to manufacture dust, 
which requires collisions to occur between its constituent planetes- 
imals. However, these collisions should not occur too frequently, 
otherwise the planetesimals will destroy themselves on a time scale 
that is smaller than the stellar age. For a planetesimal disc to sur- 
vive for an age t*, the collisional time between planetesimals needs 
to be icoii > i* . Such a condition has two regimes demarcated by 
whether the planetesimals are self-gravitating. In the regime where 



gravitational focusing becomes non-negligible, the Safronov num- 
ber O equals or exceeds unity. 



r> 



2-kGp 



1/2 



(7) 



When equation ([7| is fulfilled, we apply the following constraint. 



<Tr > 



8f2GMii,,Ur 



1/2 



{ a^ 



1/4 



(8) 



where /2 ~ 1.521 (see Appendix A of jHeng & Tremaine|2010| . 
Otherwise, we apply the constraint. 



r > 



12/iMdisc^* / GA/, 



1/2 



(9) 



where /i ^ 0.690. 

The last condition for the survival of a dynamically hot disc 
over a time t* involves gravitational scattering, which generally 
produces substantial changes in eo and iq. We demand that igrav > 
t*, which yields 

f pMdlsct ^ ^'''^ 



a,- > 



V 3/„ 



\ G^/8„-7/8 ^1/8^3/4 ^^^^,) 



/\l/4 



(10) 



The quantity Si = Si(io/eo) generally depends on the ratio of 
the rms inclination to eccentricity of the planetesimals; we have 
Si (0.5) ^ 4.50. Associated with the quantity C' = C + S2/S1, 
we have 



2C = In (A^ + 1) - In (A^ + 1) 



1 



A = 



87rGpr3 



GM-, 



1/2 



+ 



A2 

/ 87rpr" 



Ai + l' 



3\ 1/3' 



V 9A/* 



Ac = 



15cr; 
47rGpr2 



1 + 



47rpr a 
3M* 



5afa 1 fS-Kpr' 



3\ 2/3' 



+ 



2GA'/* 2 V 9M* 



In 1/2 



(11) 



Furthermore, the quantity S2 contains exponential integrals: 

S2^ l0.13[W{x)-W{'ix)], 
_ / Sttp \ ^''^ GM^ 

W (x) = exp (x) El (x) , 



(12) 



£i(a 



exp {—xx') 



dx' . 



3.1.3 Dynamically Warm Discs 

In dynamically warm discs, we require that collisions between 
the planetesimals leave the mass distribution unchanged. The latter 
condition is fulfilled when the planetesimals do not form gravita- 
tionally bound pairs (i.e., the Safronov number should be less than 
unity), 



r < 



1/2 



.SvrGp, 

The collisions are assumed to be frequent (tcou < i*), 

1/2 



r < 






(13) 



(14) 
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40 60 

wavelength (ju.m) 

Figure 2. Examples of Planck or blackbody functions with fx = l and 
f\ = 3.5, peaking at A = 22 and 77 fim, respectively. 



but non-erosive in nature, 

7 



< TV 



a 



1/4 



fmpr-Qh 



1/2 



(15) 



The binding energy per unit mass of the planetesimals is given by 



+ 



(16) 



where Qo ~ 10*-10'^ erg g"\ ro = 2 x 10* cm, a = -0.4 and 
b = 1.3 ( |Benz & Asphaug|1999[[Stewart & Leinhardt|2009"t . We 
will explore the variation of the normalization Qo as this carries an 
uncertainty of several orders of magnitude. 

An additional (and weaker) constraint is that the collisions be- 
tween the planetesimals do not result in significant viscous spread- 
ing of the disc. 



fmpr- 
6/i/6i*Mdis, 



1/2 



{GM, 



fl^a''l\ (17) 



where we have taken the typical energy loss per unit mass in a 
collision to be j^^a^. We adopt ja — \. 



3.2 Infrared Emission 

3.2.1 From Dust Grains 

Consider a dust grain, with a radius r, to be located at a dis- 
tance a from a star. It is heated to a temperature of 

T = TA'-^\ (^^^^ , (18) 



where Qabs and Qcmit are the absorption and emission efficien- 
cies of the grain, respectively. The preceding equation requires that 
the spherical grain intersects the incident starlight with a projected 
cross section of yrr^ and redistrib utes the heat over a surface area 

. |2012|, we define 



of 47rr . Following 



Lestrade et al. 



/t 



ti^abs 



1/4 



(19) 



Blackbody grains have /t = 1. Small grains may possess tempera- 
tures higher than predicted by their blackbody values (/t > 1), be- 
cause they emit inefficiently at wavelengths larger than their sizes. 
For the debris disc around GJ 581, [Lestrade et al.| (|2012|) infer 



/t = 3.5tio- ll is unknown if this is a representative value for 
the dust populations of debris discs around M stars in general. We 
note that for the debris disc around the G star 6 1 Vir, |Wyatt et al.| 
pOT2| find /t = 1.9. 

When grains are re-emitting starlight at a wavelength A, an 
effect to consider is that A may not necessarily be the peak wave- 
length of the emission profile. This can be mimicked using the fol- 
lowing expression, 

r = ^, (20) 

./A A 

where Cwion = 2897.7685 /im K is Wien's displacement constant 
and /a is a dimensionless constant. Grains with /a 7^ 1 emit at a 
peak wavelength of /a A. In the example shown in Figure|2] a grain 
with /a = 1 has a temperature of about 132 K and emits at a peak 
wavelength of 22 /im. Also shown is an example with /a = 3.5. 
In this case, the grain behaves like a blackbody emitting at a peak 
wavelength of /a A ~ 77 /im. Such an approximation allows us 
to model debris discs that emit at, e.g., A = 22 /im, but do not 
possess spectral energy distributions that peak at this wavelength. 
Note that we are not implying that [Lestrade et alT| ( [2012[ > inferred 
/a = 3.5, but rather we have picked /a = 3.5 in anticipation of 
the degeneracy between /t and /a we will next discuss. 

Given the observed wavelength A as well as the grain and stel- 
lar properties, the distance of the grain (and hence that of its parent 
planetesimals) from the star can be evaluated. 



where 



2 



/du 



T-k /dust A 



/t/a. 



(21) 



(22) 



The expression for /dust reflects a number of degeneracies inher- 
ent in this parametrization. For example, a larger value of a may be 
caused by blackbody (/t = 1) grains residing at larger distances 
emitting at a peak wavelength of /a A. It may also be caused by 
"hot" (/t > 1) grains with a peak wavelength of emission coin- 
cident with the observed wavelength (/a = 1). The "hotness" of 
a grain (/t > 1) may be offset by it emitting at a shorter peak 
wavelength (/a = 1//t), causing it to resemble a blackbody grain 
emitting at a peak wavelength of A. Strictly speaking, the use of the 
free parameters /t and /a are inconsistent with the derivations of 
equation ifTSl and Wien's law, which assume a blackbody. In the 
absence of broader knowledge on the values of these parameters, 
we consider both /dust = 1 and 3.5 in our calculations. A value of 
/dust 7^ 1 may be interpreted in any of the ways just discussed. 

Although our planetesimals are assumed to be mono-disperse 
(i.e., of a single radius/size), we allow collisions between them to 
produce a coUisional cascade of dust grains. As described in §7.4 
of [Heng & Trem aine (2010), our calculation of the infrared excess 
does not require the specification of the values of the minimum and 
maximum radii in the cascade. (Physically, a minimum radius is 
not needed because small grains do not contribute significantly to 
the infrared emission if r <C A/27r.) The infrared excess at a given 
wavelength is 



,/ir(A) = 



Bx{\T) 



K 



B,{X,n)(A-q){q-3)Rl 



(23) 



Unlike for a, there is a subtlety concerning the specification of 
the grain temperature T. It is Cwicn/T/A if one wishes to model 
hot (non-blackbody) grains emittingly predominantly at the ob- 
served wavelength A. However, it is Cwien/f\X if one wishes to 
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Figure 3. Planetesimal radius versus disc mass for a suite of model discs 
that satisfy the dynamical survival conditions stated in the text. We examine 
dynamically hot discs that are capable of generating dust via collisions. 
For illustration, we assume ti, = 300 Myr and A = 22 fim. Top panel: 
M4 star (a Ri 0.3 AU at A = 22 ^m). Bottom panel: G2 star (a ^ 4.5 
AU at A = 22 /im). Note that the labels apply only for the boundaries 
corresponding to the /dust = 1 case. 



model blackbody grains with lower temperatures being observed at 
a wavelength A shorter than the peak wavelength of emission (/a A). 
Such considerations have no effect on elucidating the allowed re- 
gion of parameter space for the planetesimals. The boundaries of 
the allowed region of /ir versus Mdisc are largely unchanged with 
the maximum value of the infrared excess, of the entire ensemble, 
differing by a factor ~ 10 between /t = 3.5 and fx = 3.5. In- 
stead, we find that the main effect of /dust is in changing the value 
of o oc /dust- Thus, to keep our analysis (and figures) simple, we 
use Cwicn/A as the input grain temperature in equation l |23| l while 
allowing for /dust 7^ 1 in equation \2l\ . 

The quantity K is the normalization factor in the radius distri- 
bution of the dust grains. 



K = 



1/2 



.iWdisc 



27r/2 Gpr^ 
3/1 0-2 



1/2 



2Q&,, 



5/12 



(24) 

where Q^ dust ~ 10^ erg g~^ is the binding energy per unit mass 
of dust grains in the collisional cascade and is approximately con- 
stant with mass (unlike for the planetesimals). The Planck function 
is given by B\ = (2/ic^/A^)/[exp(/ic/AfcBT) - 1] where h is the 
Planck constant, c is the speed of light and /cb is the Boltzmann 
constant. Denoting the number of dust grains by A'^dust, the size 
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Figure 4. Same as Figure |3] but for the infrared excess at A = 22 /^m 
associated with dust produced by colliding planetesimals in dynamically 
hot discs. For illustration, the horizontal dotted line shows the detection 
limit of /iR,thros = 0.1 associated with the WISE satellite. Note that the 
labels apply only for the boundaries corresponding to the /dust = 1 case. 



distribution of the dust grains takes the form dNdust/dr — Kr "^ 
with q = 7/2. 

3.2.2 From Planetesimals 

The planetesimals are large enough that it is safe to assume 
r/A ^ 1 and /t = 1. However, one still needs to consider the 
possibility that they may emit at a peak wavelength of f\X where 
/a 7^ 1. Thus, we set /dust = /a when using equation pTl. The 
infrared excess at a given wavelength is 



/iR = 



3Mdisc 5a (A, T) 
AnprRl Bx{X,n)' 



(25) 



Since /dust = /a, the input planetesimal temperature T takes the 
form of Cwion/,/dustA. 



4 RESULTS 

4.1 Basic Model 

Only three out of the 9 parameters in our model are either un- 
constrained by the observations or not easily specified from first 
principles: the disc mass Mdisc, the planetesimal radius r and the 
radial velocity dispersion of the planetesimal ar. We use a Monte 
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Carlo method to explore the variation of these parameters over a 
wide range: 10"^° ^ Mdisc/M® sC 1, 1 cm sC r ^ 10* km 
and 1 cm s^^ ^ a,. ^ 100 km s~^. The very wide range of pa- 
rameter values probed ensures that the allowed parameter space is 
constrained by physics rather than artificially selected boundaries. 
We assume that these parameters are uniformly distributed over the 
stated, logarithmic ranges of values. For each set of stellar param- 
eters, we generate 30000 model planetesimal discs. For each disc, 
we check the list of survival conditions described in §3.1[ if all 
of the relevant conditions are fulfilled, we deem the disk to have 
"survived". Its existence — but not necessarily its detection — is thus 
plausible. 

Figure[3]shows the suite of dynamically hot planetesimal discs 
that satisfy the dynamical survival conditions in the parameter 
space of r versus Mdisc- For illustration, we pick t, — 300 Myr and 
set the value of a to correspond to A = 22 /im as stated in equation 
pTl. We examine conditions both around a M4 (top panel) and a G2 
(bottom panel) star. As shown in Figure [3] the allowed parameter 
space is severely constrained by the survival conditions. In particu- 
lar, the icoii > i* constraint sets a lower bound on r at a given a,,, 
while demanding that the discs are dynamically hot (cr^ > (Jhot) 
and contain more than one planetesimal (A'^ > 1) set an upper 
bound. Collectively, these conditions set an upper limit on the max- 
imum planetesimal radius allowed, which is ~ 100 km for our 300 
Myr-old M4 star if /dust = 1 and the disc mass is ~ lO"^M0.The 



other survival conditions do not dictate the bounds of the allowed 
region, but rather the density of points within the region. 

One may ask if the differences in the allowed parameter re- 
gions seen for the M4 and G2 stars are due to the differences in the 
Af*, i?* and T* values or because of the difference in the value of 
a at a stated wavelength (equation ]21|). To resolve this question, 
we fix a and examine plots (not shown) similar to those in Figure 
[3] for both the M4 and G2 stars. We find that the allowed regions 
of parameter space are both qualitatively and quantitatively similar. 
Thus, the most important effect of varying the stellar type is in the 
determination of a at a fixed A, which then affects all of the pre- 
viously stated survival conditions involving a. In FigurelS] the fact 
that more massive discs are allowed around a G2 star at A = 22 
^m is primarily an effect of a ~ 4.5 AU (compared to a ~ 0.3 
AU around a M4 star) for /dust = 1- Dynamical processes tend to 
be more forgiving and act more slowly when one is located farther 
away from the star When /dust ~ 3.5, the corresponding values of 
a are larger and thus more massive planetesimal discs are generally 
allowed. 

At this point, it is worth emphasizing that the focus of our 
study is neither to examine planetesimal and debris discs around 
Sun-like stars, nor to perform detailed comparisons of synthetic in- 
frared excesses with detected ones from Sun-like stars. Applied to 
G2 stars, our models seek to demonstrate broad consistency while 
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acknowledging the degeneracies involving /dust- In other words, 
we simply wish to show that it is not difficult to produce synthetic 
infrared excesses around G2 stars that are of the same order of mag- 
nitude as in observed systems. 

Figure PI shows the infrared excess /ir, associated with dust 
grains produced in collisions between planetesimals residing in dy- 
namically hot discs. It is important to note that we are not intending 
for our calculations to be a rigorous prediction of detection statis- 
tics. A few aspects of these results are worth emphasizing. Con- 
sider the /dust = 1 case. Firstly, /ir generally increases with the 
disc mass Mdisc in a monotonic fashion. Secondly, at a fixed value 
of Mdisc, there is a dispersion in the value of /ir over several or- 
ders of magnitude. This degeneracy between the infrared excess 
and disc mass arises from the fact that for a given value of Afdisc, 
there exists a large number of solutions of r and Or that satisfy the 
dynamical survival conditions. Thus, /ir is a poor diagnostic of 
the disc mass, as one expects. Thirdly, the fcoii > t* condition sets 
an upper limit on the allowed value of /ir, because model discs 
with shorter collisions times associated with the planetesimals are 
not expected to survive or persist for the stellar age. For the 300 
Myr-old M4 star, we have /ir < 0.01, consistent with the detec- 
tion threshold stated in |Avenhaus, Schmid & MeyerH2012| for the 
WISE satellite. By contrast, we have /ir < 10 for the G2 star, 
which is broadly consistent with the measured values of /ir ~ 1- 
10 associated with Sun-like stars fBryde n et al.|2006l|Chen et al.| 
|2006[[Meyer et al.|2008^ . Thus, our survival models simultaneously 
explain the null results of Avenhaus, Schmid & Meyer (2012) in- 
volving M stars and are broadly consistent with the detections of 
previous studies examining Sun-like stars without any need for the 
finetuning of parameters. Fourthly, we note that our suite of mod- 
els suggest the existence of a vast population of planetesimal discs 
which are currently unobservable because the associated infrared 
emission is too faint to detect. For comparison, the a « 1 AU 
planetesimal disc in HD 69830 is estimated to have a disc mass 
Mdisc J; 3-4 X 10"'^ M®, several times more massive than our 
own asteroid belt (iHeng]20TTJ. 

When /dust = 3.5, the dust grains produced are either hotter 
than predicted by their blackbody temperatures at a given distance 
or emitting at a longer peak wavelength than the observed one. Both 
interpretations require the grains to be situated farther away from 
the star by a factor of /lust = 12.25. Dynamical conditions be- 
come more forgiving and the region of allowed parameter space for 
the planetesimals is larger, as seen in Figure [3] Corresponding, the 
infrared excesses allowed are typically higher as reflected in Fig- 
ure p] because the allowed disc masses are generally higher. If the 
sample of M stars examined by jAvenhaus, Schmid & MeyerH2012) 
have ages ~ 300 Myr, then a plausible, alternative interpretation is 
that they are somehow scrutinizing the fainter, undetectable mem- 
bers of the debris disc population with /dust = 3.5. 

Figure [5] shows the allowed planetesimal radius versus disc 
mass for a suite of dynamically warm discs. The condition tcoii < 
t^, plays the opposite role of tcoii > i*: it sets an upper limit on 
the planetesimal radius for a given disc mass. The maximum value 
of the planetesimal radius is set by a combination of demanding 
that the planetesimals are non-self-gravitating (0 < 1) and that 
the collisions between them are non-erosive, since both conditions 
depend on ay. The non-erosive condition restricts the rest of the 
allowed parameter space. The corresponding infrared excess from 
the planetesimals is shown in Figure [6] The upper limit to /ir for 
a wide range of A/disc is set by the minimum value of r assumed 
in our Monte Carlo calculations, but this ultimately does not affect 
the maximum value of /ir for a given suite of model discs. Rather, 



it is set by the non-erosive condition, which depends on the value 
of Qq, the binding energy of the planetesimals. The most uncer- 
tain quantity in the calculation of max{/iR} is the normalization 
Qo in equation l |16^ for the binding energy — the computed infrared 
excess scales roughly linearly with Qq. Higher values of Qo allow 
more planetesimals of a given size to be contained within a dy- 
namically warm disc and still avoid dust-manufacturing collisions, 
thereby allowing for higher disc masses and infrared excesses. 

Performing further calculations with t* = 300 Myr (not 
shown), we estimate that the predicted values of /ir are consis- 
tent with the null results of |Avenhaus, Schmid & Meyer| ^2012| l if 
Qo ^ 10^ erg g~^ and /dust = 1. For older dynamically warm 
discs, the range of Qo values allowed is less restrictive. Generally, 
our model warm discs are also broadly consistent with the reported 
values of /ir for Sun-like stars even allowing for Qo ~ 10^-10^ 
erg g~^ . Overall, if an observed disc system is a dynamically warm 
planetesimal disc, then it tends to be more massive than "tradi- 
tional " debris discs (arising from dynamically hot discs) to com- 
pensate for the relative faintness of the infrared emission from the 
planetesimals. 



4.2 Maximum Infrared Excess and Comparison to 
Observations 

The maximum infrared excess is computed across all M spec- 
tral types in FigurelTlfor tt, — 30, 100 and 300 Myr as well as 1 
Gyr. Dynamically warm discs produce similar values of max{/iR} 
to dynamically hot discs if we assume Qo = 10^ erg g"'^. As ex- 
pected, the maximum value of /ir declines with both stellar age 
and luminosity, the latter arising from the effect of probing smaller 
values of a as previously discussed. 

For comparison, we include the detection of an infrared excess 
from AU Mic at A = 22 ^m by' Avenhaus, Schmid & Meyer|p012| 
(see Spl. The data point is easily consistent with the maximum in- 
frared excess allowed for its age (i,t ~ 12 Myr), whether we have 
/dust ~ 1 or 3.5. (The "birth ring" of AU Mic is observed to reside 
at about 35 AU rather than 1 AU, cf. |Wilner et al.|2012| ) 

Figure [8] shows the maximum infrared excess across a range 
of stellar ages, t* = 0.01-10 Gyr. Generally, we have /ir oc 47^ in 
agreement with[ Dominik & Decin ( 2003 ), Wyatt et al. ( 2007 1 and 
|Heng & Tremaine, ( ,2010) for dynamically hot discs. Both dynami- 
cally hot and warm planetesimal discs around Sun-like stars gener- 
ally produce brighter infrared emission due to the more favourable 
dynamical conditions present for the planetesimals to survive for 
long time scales. We include in Figurelslthe reported values of /ir 
(at 70 pm) and f* for G stars from |Bry den et al.| l |2006^ . We do 
not specialize to the individual stellar parameters associated with 
these four objects, because doing so will only introduce minor cor- 
rections and does not change our overall result and conclusion. Our 
computed upper limits on /ir for G stars are easily consistent, by 
several orders of magnitude, with these measured values, imply- 
ing the existence of model disc solutions capable of reproducing 
the measured /ir values. For i* > 300 Myr, our models predict 
that debris discs around M stars are at best marginally detectable 
(/ir < 0.1) for spectral type M4 and later if /dust = 1. At smaller 
wavelengths (3.4 fim, 4.6 ^m and 12 /xm; not shown), our model 
discs produce lower values of the maximum allowed infrared ex- 
cess. 

We note that the /ir 3^ 1 values are implausible because 
the corresponding covering fraction associated with either the dust 
grains (for debris discs) or planetesimals (for dynamically warm 
discs) becomes optically thick. 
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Figure 7. Maximum infrared excesses (triangles) as computed by a Monte Carlo implementation of the survival models across the M spectral types. Shown is 
the data point for AU Mic (circle; t^ Ri 12 Myr), which is consistent with the maximum infrared excess allowed. Top row: dynamically hot discs, which are 
the progenitors of debris discs in the traditional sense. Bottom row: dynamically warm discs (assuming Qq = IC* erg g~^). Left column: /dust = 1- Right 
column: /dust = 3.5. 



Additionally, we compute the theoretical upper limits on /ir 
associated with AU Mic. It is apparent that these upper limits are 
again easily consistent with the detected infrared excesses at 22 
/xm and 70 /im (see SJ2l. We do not invert the measured infrared 
excesses from AU Mic to obtain an estimate for the disc mass, be- 
cause it has been previously shown that at ~ 10 Myr the disc mass 
is essentially unconstrained (|Heng|20TT|l. 



Finally, we consider the recent discovery of an old debris 
disc around the M3 star in GJ 581 via the Herschel Space Ob- 
servatory, using the PACS instrument, at 70, 100 and 160 ^m by 
|Lestrade et al.|pbl2 i. We use the stellar parameters listed in Ta- 
ble 2 of ' Lestrade et al.| ( |2012| l and take the average of the range 
of stellar ages quoted (5 Gyr). At 70, 100 and 160 /im, we ob- 
tain max{/iR} « 3, 25 and 200, respectively, for dynamically hot 
discs if /dust = 1- The values of max{/iR} become even larger for 
/dust > 1. For dynamically warm discs, we obtain max{/iR,} ~ 9, 
80 and 900 at the same respective wavelengths for Qq — 10^ erg 
g~^ and /dust = 1; these upper limits are higher for larger values 
of Qo and /dust. These computed values of the maximum infrared 
excesses are consistent with the values of about 3, 7 and 20 reported 
in|Lestrade et al.|((2012|). (See their Figure 6.) 



5 DISCUSSION 

5.1 Summary 

We have used survival models to examine the occurrence of 
planetesimal and debris discs around M and G stars. The salient 
points of our investigation include: 

• The dynamical survival conditions severely restrict the range 
of disc masses and planetesimal sizes allowed for a disc to persist 
for its stellar age. Despite this, the infrared excess is a poor diagnos- 
tic of the disc mass at young ages (~ 300 Myr), generally spanning 
several orders of magnitude for a given disc massrl 

• The dearth of detectable infrared excesses from M stars, as 
probed by the WISE satellite, is possibly due to the small semi- 
major axes probed at 3.4-22 fim. These small semi-major axes 
translate to less forgiving dynamical survival conditions for the 
planetesimal discs. In other words, planetesimals may not exist at 
small distances from the star for time scales > 300 Myr, regardless 
of whether the spectral type is M or G. This interpretation becomes 
less clear when the dust grains (in debris discs) or planetesimals (in 
dynamically warm discs) deviate from blackbody behaviour (i.e., 
are hotter) or emit predominantly at wavelengths longer than the 



■^ But see |Heng|J201l| for an application of the survival models to estimate 
the disc mass when t* ~ 1 Gyr. 
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observed one. However, the interpretation is acceptable when they 
are hotter than blackbody and emit at peak wavelengths shorter 
than the observed one. Clearly, some of these degeneracies may be 
broken by searching for debris discs around M stars across a broad 
range of wavelengths. 

• Both dynamically hot and warm planetesimal discs are capa- 
ble of producing infrared excesses that are consistent with the re- 
ported values for G stars. It is plausible that both dynamically hot 
and warm planetesimal discs will generate detectable infrared ex- 
cess around M stars as well, the former through dust grains gener- 
ated during planetesimal collisions and the latter via the planetesi- 
mals reprocessing the starlight themselves. 

• The detected infrared excesses from AU Mic and GJ 581 are 
easily reconciled with youth and the large distances probed by the 
observations, respectively. These conclusions are unaffected by the 
degeneracies involving the non-blackbody behaviour (or peak ver- 
sus observed wavelength probed) of dust grains. 



5.2 Are Non-Gravitational Forces Important? 

A reasonable concern with our analysis is that we have ne- 
glected the effects of non-gravitational forces. For Sun-like stars, 
|Heng & Tremaine, ^2010^ previously demonstrated that Poynting- 



Robertson drag may be neglected when computing the emission 
properties of dust grains (see their §5). However, this may not gen- 
erally be the case for M stars, since |Plavchan et al.| ( |200 5) and 
|Strubbe & Chiang ] (|2006|l have shown that stellar wind drag may 
dominate over Poynting-Robertson drag. This result may be repro- 
duced by examining the efficiency factor describing both types of 
non-gravitational drag. 






(26) 



where Qpr < 1 is the efficiency factor of Poynting-Robertson 
drag. When 2-Kr/\ <C 1, we have Qpr ~ 2-nr/\. We estimate that 



the stellar luminosity for M dwarfs is Ct — 4:TTR^asBT^ 



10" 



10 ^JZq with Cq denoting the solar luminosity. The stellar mass 



loss rate for M stars is A/* ~ lO-lO^Mrr) with M, 



denoting the solar mass loss rate (see 



2L. 



Plavchan et al 



1.3xl0 ^^g 
and 



2005 



references therein). Essentially, we infer that 



M* > Mq and £* < Cc, 



M* 



> Qpr. 



(27) 



Thus, we have Qdrag ~ Mi,(? / d, for M stars, consistent with the 
conclusions of |Plavchan et aI.H2005| and |Strubbe & Chiang| ( [2006| >. 
That stellar wind drag acts much faster than Poynting- 
Robertson drag, for M stars, does not address the issue of whether 
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inter-grain collisions occur on an even faster time scale, in which 
case they are destroyed before experiencing significant orbital de- 
cay (iStmbbe & Chiang 2006 ). For Sun-lilse stars,'Wyatt'("2005') and 
|Heng & Tremaine (20101 showed that collisions generally domi- 
nate over Poynting-Robertson drag. We now elucidate the condi- 
tions under which collisions dominate over stellar wind drag for 
dust populations around M stars. 

Denoting the radius of a dust grain by r, the time scale for 
Poynting-Robertson and stellar-wind drag to act on it is ( [Bums et] 
|al.|1979| i 



idr 



Sttc a pr Sna pr 



3Ci,Qdi 



3Af* 



(28) 



The collisional time scale associated with the dust grains is ( |Heng| 
[&TVemaine 2010) 
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Stellar wind drag may be neglected (icoii < tdrag) when 
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An alternative way of understanding equation ((30b is to recast 
it using equation \2\\, which yields an inequality for the observed 
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(32) 



Generally, observing a debris disc at longer wavelengths corre- 
sponds to examining it at larger distances from the star. Collisions 
occur on a time scale of fcoii oc a^", while stellar wind drag occurs 
on a time scale of tdrag oc a^, implying that there is a maximum 
distance beyond which collisions do not occur fast enough. This 
can be compensated by having a more massive disc (higher Afdisc), 
a lower stellar mass loss rate (lower A/*) or a cooler star (lower 
Tt). All of these properties are reflected in equation ( |31^ . 

Yet another approach is to derive the minimum disc mass 
above which stellar wind drag may be neglected (i.e., Afdisc > 



A^di. 



.), 



Afdi: 



disc,sw — 



Tr/mAf* / Aswr*/, 



3/2 



(GAL) 



'1/2 



32/l \ Cwicn 

(33) 
Essentially, by estimating values for A/disc, sw we can evaluate the 
regions of parameter space, as shown in our figures, where stellar 
wind drag does not affect our computed results. 

Figure |9] shows calculations for Asw as a function of Afdisc 
across the stellar types MO to M9. Given the large uncertainty in 
the stellar mass loss rate for M stars, we include calculations for 
both Af* = WMq and W^Mq. It is apparent that for A/disc > 
10"*Afa), we obtain Asw ~ 100-1000 ^m for /dust = 1; for 
/dust = 3.5, Asw decreases by a factor of /dust- These estimates 
imply that for debris discs that are at least roughly as massive as 
our asteroid belt, the effects of stellar wind drag may be neglected 
when they are observed at infrared wavelengths. More precise con- 
straints may be obtained by examining the values of Afdisc, sw, also 
shown in Figure[9 To relate these estimates to both the observations 
of [Plavchan et at ( |2005[ l and |Avenhaus, Schmid & Meyer|p012| ), 



we adopt Asw = 12 and 22 /im, respectively. Generally, the min 
imum disc mass above which stellar wind drag may be neglected 



10- 



-10 



/dust-^ffi 



Plavchan et al. 



2005 I 



is low: A'/disc,sw 

searched for A = 12 pm infrared excesses around 9 M stars, but 
did not report any detections. If stellar wind drag is to be invoked 
as an explanation for these non-detections, the discs around these 
M stars must have masses lower than the A/disc, sw values estimated 
at Asw = 12 pm (even leaving aside the issue that such low-mass 
discs probably emit infrared radiation below the detection limits). 

Among the ensembles of both debris discs and dynamically 
warm planetesimal discs simulated by our Monte Carlo calcula- 
tions, the most massive members are unaffected by stellar wind 
drag. Even considering the wide range of stellar mass loss rates as- 
sociated with M stars, the discs that are affected by stellar wind 
drag have low enough masses that they are undetectable using cur- 
rent instrumentation. When stellar wind drag is considered for these 
low-mass discs, the associated dust grains (or planetesimals, if con- 
sidering dynamically warm discs) become both undetectable and 
non-existent. 



5.3 Future Prospects 

Our results in Figures 17] and [S] suggest that M stars with 
t* > 300 Myr are at best marginally detectable by the WISE 
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satellite at 22 fim, while those with i* ~ 10 Myr should be eas- 
ily detectable — an age effect. The caveat to this statement is that 
the dust grains (in debris discs) or planetesimals (in dynamically 
warm discs) either behave like or mimic blackbodies emitting pre- 
dominantly at the observed wavelength (see j ]3.2| for details). Con- 
versely, if this caveat was strongly violated, then Figures [7] and [8] 
predict that debris discs around M stars should be prevalent, a phe- 
nomenon that is not observed. Considerations of formation, rather 
than survival, may further inform this issue. 

The ability to compute /ir values that are consistent with the 
measured ones for the 12 Myr-old AU Mic debris disc further but- 
tresses the suggestion about the non-detections by WISE being an 
age effect. Gas dispersal in protoplanetary discs around M stars 
is expected to take between 8 and 12 Myr fSimon et ^2012i, 
somewhat longer than for Sun-like stars. [Schneider, Metis & Song| 
( |2012^ examined 30 K and M stars in the TW Hydrae association 
(ti, « 8 Myr) and found that 42tg°% of them have infrared ex- 
cesses of /iR ~ 0.1-10 (using the 22 ^m channel of WISE), in- 
dicating the presence of dusty protoplanetary discs. [Simon et al.| 
( [2012 ) discovered a 6 Myr-old, pre-main-sequence, M4 star with 
/iR « 40 (again using the 22 /im channel of WISE) in the rj Cha 
star cluster. (See also [Gautier et al.|2 008 ) No infrared excess was 
detected in the older star clusters: Tuc-Hor (i* ~ 30 Myr) and 
AB Dor (t* ~ 70 Myr). Taken together, these results suggest that 
young (t* ~ 10 Myr) M stars should host debris discs that are 
detectable with WISE. 



Plavchan, P., Werner, M.W., Chen, C.H., Stapelfeldt, K.R., Su, 

K.Y.L., Stauffer, J.R., & Song, I. 2009, ApJ, 698, 1068 
Rieke, G.H., et al. 2005, ApJ, 620, 1010 
Rizzuto, A.C., Ireland, M.J., & Zucker, D.B. 2012, MNRAS, 421, 

L97 
Schneider, A., Mehs, C. & Song, I. 2012, ApJ, 754, 39 
Simon, M., Schlieder, I.E., Constantin, A.-M., & Silverstein, M. 

2012, ApJ, 751, 114 
Stewart, S.T., & Leinhardt, Z.M. 2009, ApJ, 691, L133 
Strubbe, L.E., & Chiang, E.I. 2006, ApJ, 648, 652 
Su, K.Y.L., et al. 2006, ApJ, 653, 675 
Toomre, A. 1964, ApJ, 139, 1217 
Urban, L.E., Rieke, G., Su, K., & Trilling, D.E. 2012, ApJ, 750, 

98 
Wilner, D.J., Andrews, S.M., MacGregor, M.A., & Hughes, A.M. 

2012,ApJL, 749, L27 
Wyatt, M.C. 2005, A&A, 433, 1007 
Wyatt, M.C, Smith, R., Su, K.Y.L., Rieke, G.H., & Greaves, J.S. 

2007, ApJ, 663, 365 
Wyatt, M.C. 2008, ARA&A, 46, 339 
Wyatt, M.C, et al. 2012, MNRAS, 424, 1206 
Zuckerman, B. 2001, ARA&A, 39, 549 



Acknowledgments KH acknowledges a tenure-track assistant pro- 
fessorship from the Center for Space and Habitability (CSH) at 
the University of Bern, an affiliate membership at the Institute for 
Theoretical Physics (ITP) at the University of Ziirich and partial 
financial support from the Swiss-based MERAC Foundation. KH 
conducted part of the research while holding a Zwicky Prize Fel- 
lowship from ETH Ziirich in 2012. We thank Michael Meyer, Hans 
Martin Schmid and Henning Avenhaus for interesting discussions. 
We are grateful to the anonymous referee for constructive reports 
that improved the quality and clarity of the manuscript. 



REFERENCES 

Avenhaus, H., Schmid, H.M., & Meyer, M.R. 2012, A&A, 548, 

A105 
Backman, D.E., & Paresce, F. 1993, in Protostars and Planets III, 

p. 1253-1304 
Benz, W, & Asphaug, E. 1999, Icarus, 142, 5 
Bryden, G., et al. 2006, ApJ, 636, 1098 
Bums, J.A., Lamy, PL., & Soter, S. 1979, Icarus, 40, 1 
Chen, C.H., et al. 2006, ApJ, 166, 351 
Dominik, C, & Decin, G. 2003, ApJ, 598, 626 
Gautier III, T.N., RebuU, L.M., Stapelfeldt, K.R., & Mainzer, A. 

2008, ApJ, 683, 813 
Heng, K., & Tremaine, S. 2010, MNRAS, 401, 867 
Heng, K. 201 1, MNRAS, 415, 3365 
Kaltenegger, L., & Traub, W.A. 2009, ApJ, 698, 519 
Kenyon, S.J., & Bromley, B.C. 2010, ApJS, 188, 242 
Krivov, A.V., MuUer, S., Lohne, T, & Mutschke, H. 2008, ApJ, 

687, 608 
Lestrade, J.-F, et al. 2012, A&A, 548, A86 
Meyer, M.R., et al. 2008, ApJ, 673, L181 
Plavchan, P, Jura, M., & Lipscy, S.J. 2005, ApJ, 631, 1 161 



© 2013 RAS, MNRAS 000.[T|pT] 



